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Abstract. 

This work is a direct and model-independent measurement of the low mass n + n~ phase motion 
in the D + — ► n~ K + K + decay. The results show a strong phase variation, compatible with an isoscalar 
(7 (500) meson. This result confirms the previous Fermilab E791 result which found evidence for 
the existence of this scalar particle using a full Dalitz-plot analysis. 



INTRODUCTION 

Charm meson decays are a natural place for studying light scalar mesons. These decays 
are a clean environment: well defined initial state, very small non-resonant component 
and large coupling to scalars. In addition, D decays can provide insight into the quark 
content of these controversial particles, since the bulk of the hadronic width comes from 
modes for which there is a W-radiation amplitude. 

Studies of light scalars using charm decays started a few years ago, when E791[l] 
presented results of Dalitz-plot analyses of the decays Z)+,Z) + — ► n~n + n + [2, 3] and 
D + — > K~n + Tt + [4]. In particular, the analysis of the decays D + — > tz^tz + tz + and D + — > 
K~tz + tz + showed evidence for the o and k mesons. In the case of D + — > n~n + n + , the 
o appears as an accumulation of signal events in low iz + iz~ mass. The E791 data could 
only be described by an amplitude having an s dependent phase (s being the mass 
squared). 

Several years later, the same effect - an accumulation of signal events in low iz + n~ 
mass - was also observed in other decays from different experiments, for, instance, 
D° -> K q k + k- from CLEO[5] and Belle[6], and J/y -> (on + n- from BES[7]. Again, 
the description of these data requires amplitudes with s-dependent phase. 

Since E791 publication of the — > tz^tz + tz + results, two kinds of criticisms 
have been made: the quoted values of both o and K parameters are not correct because 
the simple Breit-Wigner formula is inadequate to describe broad scalars, especially when 
near the threshold; the other criticism is that one should not claim the existence of a 
resonance without showing the phase motion of the corresponding amplitude. 

The Breit-Wigner formula, although being in this case only a naive approximation, has 
the key ingredient: an s-dependent phase. On the other hand, there is no single agreed 
way to treat broad scalars near threshold. The o parameters depend strongly on the 
assumed functional form of its line shape. 

In this work[8] the second criticism is addressed. The phase variation of the low-mass 
7T + 7T~ amplitude is extracted in a model independent way. 



THE AMPLITUDE DIFFERENCE METHOD 



In Dalitz-plot analysis, resonant amplitudes are complex functions written in the general 
form srf = f(s)e l8 ( s \ Non-resonant amplitudes, in contrast, have 8(s) = constant. If 
two resonant amplitudes cross in some region of the Dalitz-plot, they will interfere. The 
interference pattern in this crossing region depends on the s-dependent phases of both 
resonances. One can extract the phase motion of an unknown amplitude that crosses a 
well known resonance provided that: 

• the contribution of other amplitudes is negligible in the crossing region between the 
amplitude under study and the known resonance; 

• the integrated amplitude of the known resonance is symmetric with respect to an 
effective mass squared (tn 2 e ff)- 

The second condition ensures that, by comparing the amplitude below and above 
m^ff, we end up with an expression that involves only the desired phase 5(513). Then, 
we can write the approximate amplitude of this phase space region in a simple way, 

^{sn,s u ) - a R ®W{s n ) ^(s n ,s 13 ) + a,/(p*Mi 3 ) sind(s 13 ) ^ 5 ^)+r) (1) 

where 7 is the overall relative final state interaction (FSI) phase difference between the 
two amplitudes, cir and a s are respectively the real magnitudes of the known resonance 
and the under-study complex amplitude, sin8(si3)e' 8 ( Sl ^ represents the most general 
amplitude for a two-body elastic scattering; p* /y/s l3 is a phase space factor to make this 
description compatible with kk scattering and ^(512,513), SSW{s\2) are the angular 
function and Breit-Wigner for the known resonance, respectively. 

ThequantityA | srf(s\i) | 2 =| srf '( m e// + £ ' 5 i3) I 2 ~~ I ^( m lff ~~ e ' s i3) | 2 , which is the 
difference of the amplitudes squared after integration over syi, is computed in bins of 
s\ 3. It takes the form 

2= -4a^/(pV^ 3 ) g mor {sin{l8{su) + y) _ siny) ^ {su y {pV ^ s } 
£2 + m o r o 

(2) 



If 5(513) is an analytical function of 513, then there will be maximum and minimum 
values of A | srf (513) | 2 , which we can use to determine both the constant term in the 
above equation and the phase 7, 

- wo. w ye-or. sV = (A , ^ lL> _ A , lL)/2 (3) 

e + m o r o 

y=sin ( a h 2 -a h4 r ) (4) 

^ I \min ^ I *** Imax 

where srf' = srf j{Jt 'y/s n / 'p*). 
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FIGURE 1. D + — > n~n + n + Dalitz-plot folded distribution. The events used by the AD analysis are in 
the hatched region. The size of the area of each bin in the plot corresponds to the number of events in that 
bin. 



Finally, considering 5(513) an increasing function of 513, we have 



This is, in essence, the idea of the amplitude difference method [9]. The method 
was shown to work in a "calibration" exercise using the /o(980) resonance in the 
Df — > n~n + 7t + decay [10]. In this case we have the /o(980) contribution in both 512 
and *i3 axes. We were able to get a phase motion 5(513) compatible with the /o(980) 
using the /o(980) in 512. 



PHASE MOTION OF SCALAR LOW-MASS %~%+ AMPLITUDE 



The folded Dalitz-plot distribution of the D + — > n~n + K + decay is shown Fig. 2. The 
horizontal and vertical axes are the squares of the k + k~ invariant mass high (512) and 
low (513) combinations. 

To study the low mass region in 513, there are three possible well known resonances 
in 512 to act as a probe in this decay: p(770), /o(980) and /2(1270). Figure 1 shows 
that the p(770) and /o(980) are located in regions where other amplitudes can not be 
considered negligible. On the other hand, the tensor /2(1270), m 2 , = 1.61 GeV 2 /c 4 , is 
placed where the p(770), in the crossed channel reaches a minimum due to its decay 
angular distribution. 

With the proper choice of i%ff = 1-535 GeV 2 /c 4 , the integral over 512 of the fi ( 1270) 
amplitude squared is symmetrical: the number of events between m^ff and m^ff + £ (£ = 
0.26 GeV 2 /c 4 ) is equal to the number of events between m 2 . <■ and m 2 , <• — £. Moreover, 



5(5i 3 ) = ^(sin \-^A I ^(513) | 2 +sin(r))-r) 



(5) 



in this region there is no significant contribution other than the KK complex amplitude 
under study in 5 13 (the amount of p(770) within this mass region was estimated to be 
~5%). The choice of the /2(1270) as the analyser amplitude satisfies the necessary 
conditions for the amplitude difference method. 

The acceptance and the background must be similar between irf,ff and rn 2 e ff + £ and 

m^ff and rn^ff — £, otherwise there would be biases in 6(513) . Monte Carlo simulations 
show that the acceptance is nearly uniform in this region. The background in this re- 
gion comes mostly from random combinations of three pions, and it is also uniformly 
distributed. Since we are subtracting two similar distributions, we considered the back- 
ground only in the size of the statistical error. 

The /2(1270) angular function has a zero at about 513 ~0.48 GeV 2 /c 4 . This means a 
singularity in srf' . This singularity is handled in the following way. The data is divided 
into ten 513 bins. The binning is such that the singularity is placed in the middle of one 
bin. Doing this, we isolate the singularity in a single bin (bin 6) and discard its further 
use in the analysis. 

The values of 7 and ^ were obtained solving Equations 3 and 4. The value of the phase 
7 from the amplitude difference method is in agreement with that of the full Dalitz-plot 
analysis, y AD = 2.78 ±0.38 ±0.40 and jDaiitz = 2.59±0.19[3]. 

Having 7 and % ', the phase 5(513) is obtained for each 513 bin using Equation 5. 
There are ambiguities arising from the sin~ l operation, so 5(513) was determined with 
the assumption that the phase difference starts at zero at threshold and is an increasing, 
monotonic, smooth function of 513. 

The phase motion of the low mass amplitude, including systematic and statisti- 

cal errors, is shown in Fig. 2. In spite of the limited statistics, a strong phase variation is 
clearly observed. Starting from zero at the threshold the phase varies by about 180° and 
saturates at around 513 = 0.6 GeV 2 /c 4 . This is the expected behavior of resonance. The 
observed phase motion supports the interpretation of the a(500) as a true resonance. A 
constant or slowly varying phase would disfavor this interpretation. With more statis- 
tics we could have more bins and the pole position could be inferred. Fig. 2 shows also 
the phase motion of the simple Breit-Wigner (solid line) used in [3]. Even considering 
the Breit-Wigner as a naive approximation, there is a qualitative agreement between its 
phase motion and the directly extracted 5 (513): the description of the data requires an 
amplitude with a strong phase variation. 



CONCLUSIONS 

A direct, model-independent measurement of the phase motion of the low mass n + n~ 
scalar amplitude was discussed. Using the well known /2(1270) tensor meson in the 
crossing channel as the base resonance, from the D + — > k~k + k + decay, the 5(513) 
phase motion was extracted. We obtain a 5(513) variation of about 180°, which is the 
expected behavior of a resonant amplitude. This result supports the interpretation of 
the a(500) as a true resonance, in agreement with the conclusions from our previous 
analysis of full D + — > k~k + k + Dalitz-plot. We could not extract the a pole position 
with this method due to the limited statistics. The measurement of the correct o pole 




FIGURE 2. Phase motion vs 513, with errors shown (systematic and statistical in quadrature). The 
continuous line is the Breit-Wigner phase motion, with the E791 parameters for the a(500). 



position from a full Dalitz-plot analysis needs the correct functional form (theorists 
should agree on what it is). Whether the o pole is the same in charm decay and scattering 
remains an open question. 
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